Anisotropic Sliding Dynamics, Peak Effect, and Metastability in Stripe Systems 
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A variety of soft and hard condensed matter systems are known to form stripe patterns. Here we 
use numerical simulations to analyze how such stripe states depin and slide when interacting with 
a random substrate and with driving in different directions with respect to the orientation of the 
stripes. Depending on the strength and density of the substrate disorder, we find that there can be 
pronounced anisotropy in the transport produced by different dynamical flow phases. We also find 
a disorder-induced "peak effect" similar to that observed for superconducting vortex systems, which 
is marked by a transition from elastic depinning to a state where the stripe structure fragments 
or partially disorders at depinning. Under the sudden application of a driving force, we observe 
pronounced metastability effects similar to those found near the order-disorder transition associated 
with the peak effect regime for three-dimensional superconducting vortices. The characteristic 
transient time required for the system to reach a steady state diverges in the region where the flow 
changes from elastic to disordered. We also find that anisotropy of the flow persists in the presence 
of thermal disorder when thermally-induced particle hopping along the stripes dominates. The 
thermal effects can wash out the effects of the quenched disorder, leading to a thermally-induced 
stripe state. We map out the dynamical phase diagram for this system, and discuss how our results 
could be explored in electron liquid crystal systems, type-1.5 superconductors, and pattern-forming 
colloidal assemblies. 

PACS numbers: 64.60.Cn,61.20.Gy,05.40.-a 
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I. INTRODUCTION 

Stripe formation occurs in a wide variety of softi~— 
and harcUi~— condensed matter systems. These stripe 
patterns are often a consequence of some form of effec- 
tive competing or multiple length scales in the pairwisc 
interactions between the particles^ 3 - - — i 9 ' 10 i 14 ~ 16 ' 18 . For 
soft condensed matter, pattern formation can occur when 
the particles experience intermediate range repulsion and 
short range attraction, such as in certain types of col- 
loidal systems^ 9 -. In addition to stripe phases, numerous 
other patterns can appear as a function of density, tem- 
perature, or particle interaction strength, including bub- 
ble, clump, and uniform crystalline phases^ 3 - - — i 10 ' 14 : 17 : 19 . 
The competing interactions that produce the stripes may 
be produced by particles that have both a short range 
attraction and a long range repulsio n 10 i 14 i 18 i 19 ; however, 
systems with only repulsive interactions can also exhibit 
stripe phases^ 5 -^ provided that there are at least two 
length scales in the interaction potential. Typically, as 
the density increases, the system progresses from a low 
density clump phase to an intermediate density stripe 
phase, and then to a higher density bubble phase where 
organized voids appear in the system; finally, at the 
highest densities, the particles form a uniform crystal 
state^&Sii 9 -. In two-dimensional (2D) systems of finite 
size, a stripe phase containing oriented stripes is often 
observe d 4 ! 5 ! 19 ; however, for larger systems, the strong 
degeneracy in the stripe ground state orientation can 
produce a labyrinth pattern composed of many differ- 
ent stripe orientations^! 8 -^. The presence of any type 
of bias produced by the boundaries, a substrate, or an 



external drive such as a shear breaks the symmetry of 
the stripe ground state and causes the stripes to align in 
a single direction 9 - ! 20 ! 21 . 

In addition to soft matter systems, there is growing 
evidence that stripe and bubble phases occur in hard 
condensed matter systems such as 2D electrons in the 
quantum Hall regim o 17 ! 22 ' 23 and charge ordering in high 
temperature superconductors^— . Evidence for stripe 
phases in two-dimensional electron gas (2DEG) systems 
includes anisotropic transport curves which have been 
interpreted as indicating that the stripes have a single 
preferred orientation 2 ^—. If the stripes take the form 
of a charge ordered state, the transport anisotropy im- 
plies that the stripes can slide more easily when the 
drive is applied parallel to the stripes than when it is 
applied perpendicular to them. The alignment of the 
stripes in the 2DEG systems may be due to small intrin- 
sic biases that form during sample growth**. There have 
also been recent 2DEG experiments that show that dc 
drives or other external driving can dynamically orient 
the stripes under certain condition a 27 ! 28 . Other recent 
experiments have shown that the stripe direction can 
be controlled with a strain, making it possible to alter 
the anisotropy with a strain field^ 9 .. Transport experi- 
ments in 2DEGs have revealed sharp conduction thresh- 
olds, a series of intricate jumps in the current versus re- 
sistance curves, pronounced hysteresis, and changes in 
the conduction noise, suggesting that these systems are 
undergoing depinning transitions and dynamic changes 
in the sliding dynamic o 27 ' 30 ! 31 . Additional evidence for 
charge ordered states in 2DEGs has come from resonance 
measurements 32,33 . 
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Another recently described system where stripe pat- 
terns occur is in "type-1.5" superconductors, predicted to 
appear in two-band superconductors such as MgB^. In 
a type-II superconductor under a magnetic field, the flux 
in the sample takes the form of quantized vortices which 
organize into a uniform triangular lattice as a result of 
their repulsive interactions. In contrast, in type-1.5 su- 
perconductors, the vortices have both an attractive and 
a repulsive component to their interaction s 34 ' 35 , which 
in principle will lead to the formation of clumps and 
stripes. Experiments in the two-band superconductor 
systems have revealed evidence for disordered clump-like 
vortex structures; however, strong pinning in the sam- 
ples probably prevents the detection of ordered patterned 
structures 3 ^. 

In the 2D electron systems and the type-1.5 super- 
conductors, the interplay between the disorder in the 
sample and an external drive should produce very rich 
dynamics with different types of depinning transitions 
and sliding states. It should also be possible to subject 
stripe-forming soft matter systems to both an external 
drive and quenched disorder. Experiments have already 
been conducted on the depinning of purely repulsive col- 
loids interacting with quenched disorder—. Similar ex- 
periments could be performed with colloidal systems that 
have interactions which lead to stripe formation. The 
work we describe here is also relevant to systems ex- 
hibiting anisotropic sliding friction due to the formation 
of stripe-like surface ordering 3 ^. To address how stripe- 
forming systems behave in the presence of both driving 
and quenched disorder, we simulate a collection of parti- 
cles driven over randomly placed attractive pinning sites 
and interacting with a long range repulsion and a short 
range attractio n 14 ' 19 ' 21 . We show that when the stripes 
have a specific orientation, a number of distinct sliding 
states can occur which have distinct anisotropic transport 
signatures, including a peak in the anisotropy produced 
when the dynamics is plastic for driving transverse to the 
stripes but elastic for driving along the stripes. 

In our previous work, we have examined the depin- 
ning and sliding of disordered bubble, clump, and stripe 
phases, and found that for a fixed pinning density, the 
stripe phase has the highest depinning threshol d 14 ' 21 . In 
this case, the stripes did not have a single orientation 
but instead formed a disordered labyrinth pattern. When 
plastic depinning occurred, it was possible to induce a dy- 
namical reordering transition into a stripe state aligned 
with the driving directional. We also found that the dy- 
namically induced reorientation strongly depends on the 
strength of the quenched disorder. Only for sufficiently 
strong quenched disorder are there enough plastic distor- 
tions to permit the formation of the oriented stripes! 4 -. 
When the pinning is weak, the labyrinth structures depin 
elastically without any distortions and the aligned stripes 
never form. 

Here we analyze the transition from elastic to plastic 
depinning and show that for some parameters, the stripe 
system exhibits a peak effect phenomenon similar to that 



observed at the transition from elastic to plastic depin- 
ning in vortex matter. The vortex peak effect is associ- 
ated with a sharp increase in the depinning force as well 
as changes in the transport curves 3 ^—. We show that 
the peak effect in the stripe system can occur for driving 
in either direction and that it is possible to have a peak 
effect for one direction of drive but not the other. In pre- 
vious work, we showed that there is a broad maximum 
in the depinning force for the stripe phase as function 
of the strength of the attractive term. In this work, we 
study the peak effect in the stripe phase as a function 
of disorder and temperature, and find that it occurs as 
a sharp, first-order-like transition which is similar to the 
peak effect observed in superconductors. 

In this work we explicitly focus on the case where the 
stripes are already in an aligned state rather than in a 
disordered labyrinth phase. This permits us to apply a 
drive in two well-defined directions, along and perpen- 
dicular to the stripes, and to compare the anisotropic 
response for different strengths of quenched disorder. As 
noted previously, many of the 2DEG stripe systems ap- 
pear to contain oriented stripes. To our knowledge the 
depinning and sliding dynamics of an oriented stripe sys- 
tem has not previously been numerically studied. We 
find several new types of sliding phases that do not ap- 
pear for sliding dynamics in isotropic systems such as 
vortices 4 ^—, colloids 4 ^, sliding charge density waves^ 9 -, 
or sliding Wigner crystals 5 ^. For example, we find several 
different types of plastic stripe flow. In one state, indi- 
vidual structures slide past stationary stripes; in another 
state, the stripe structure remains intact but a portion 
of the particles within the stripes are pinned while other 
particles flow past in one-dimensional (ID) channels. For 
strong disorder the stripe structure breaks apart and the 
flow is similar to the plastic flow observed in isotropic 
vortex system s 44 ' 45 ' 48 . We also find that the extent to 
which the stripes reorient is strongly sweep rate depen- 
dent. This affects measurements of the depinning thresh- 
olds and features in the transport curves. 



II. SIMULATION 

In Fig. 1 we show a snapshot of our system containing 
stripes oriented in the y-direction. The easy driving di- 
rection is along y, parallel to the stripes, while the hard 
driving direction is along x, perpendicular to the stripe 
pattern. Our simulation box has periodic boundary con- 
ditions with sides L y = L and L x — 1.097L. We consider 
N = 380 particles with a density given by p — N/ (L x L y ). 
Here we fix p = 0.36. The particles interact with a 
long range Coulomb repulsion and a short range expo- 
nential attraction. The resulting interaction potential is 
repulsive at very short ranges due to the Coulomb term, 
attractive at intermediate range, and repulsive at long 
range. The dynamics of the particles are determined by 



3 



34 



y 












FIG. 1: The particle positions (black dots) for a system with 
competing long-range repulsion and short-range attraction. 
For this density of p = 0.36 we obtain a stripe state with the 
stripes aligned along the y-direction. The external drives are 
applied either along the stripe direction, , or perpendicular 
to the stripes, Ff). 

integrating the following equation of motion: 

^ = E WiRii) + Ff + Ff c + Ff. (1) 

Here Ri is the position of particle i. 

The first term on the right hand side of Eq. (1) is the 
particle-particle interaction potential 

V(R ij ) = ^--Bexp(- K R ij ) (2) 

with Rij = |Rj — Rj|, B = 2.0, and k = 1.0. To avoid the 
divergence from the Coulomb term at small Rij we place 
a constant- force cutoff at Rij < 0.1. The Coulomb term 
does not permit a long range interaction cutoff so for 
computational efficiency we employ a Lekner summation 
method to calculate the long range Coulomb forced. The 
second term of the interaction potential is a phenomeno- 
logical short range attractive interaction. 

The second term on the right in Eq. (1) is the 
force from the quenched disorder, modeled as N p non- 
overlapping randomly placed parabolic pinning sites with 
density p p — N p / (L x L y ) and with 

Ff = £ (F v /Rp) R^S(R P - fl«)*<« (3) 



Here, RL is the location of pinning site k, R p is the 
pinning radius which is set to R p = 0.2 unless other- 
wise noted, F p is the maximum force from a pinning site, 

B$ = |R* - R[ rt |, R^ 3 = (Ri - Ri p) )/^\ and Q is 
the Heaviside step function. 

The force Ff c in Eq. 1 arises from an external dc drive 
applied unidirectionally to all the particles in either the y 
or x direction, Ff c — F^y or Ff c = Fj^x. We measure 
the depinning threshold and transport curves for each 
driving direction by summing over the velocities of the 
particles, (V a ) — • a with a — x,y. 

The final term on the right hand side of Eq. (1) rep- 
resents the forces from randomly distributed thermal 
kicks with the following properties: (F?(t)) — and 
(Ff(t)Ff(t')) = 2r]k B Sij(t - t'), where k B is the Boltz- 
mann constant. In previous equilibrium studies of this 
system, we identified the densities at which different 
clump, stripe, and bubble phases occur~. Here we work 
at p = 0.36 corresponding to the case of stripes contain- 
ing approximately three particles per row as shown in 
Fig. 1. The initial particle positions were obtained from a 
very slow simulated annealing from a high T to T = 0.0. 
The stripes align in the y-direction during the anneal, 
and the pinning potential is not applied until after the 
annealing process. 

III. PULSE MEASUREMENTS AND 
METASTABILITY FOR STRONG DISORDER 

We first examine the dynamical response when differ- 
ent strengths of external drive are suddenly applied to 
the system. When the external drive is slowly increased 
from zero, the system passes through several different dy- 
namical phases. In contrast, for the sudden pulse drive, 
the system can pass directly from a pinned state to a 
sliding state, and an ordered moving state may appear 
that cannot be reached by slowly increasing the driving 
force. For strong pinning, the stripe structure breaks 
up or fragments near depinning, and the anisotropy of 
the two driving directions is reduced for a slow ramp of 
the driving force; however, for the pulse measurements, 
a pronounced anisotropy can be preserved. 

We conduct a series of pulse drive simulations at 
F p = 0.9 and p p = 0.38. For these parameters, slow 
driving ramps would produce a breakup of the stripe 
structure for driving in either the x or y direction. Af- 
ter we apply the pulse drive, the system typically passes 
through a transient state and the velocity relaxes to a 
steady state value after a characteristic time r. We con- 
struct a pulsed-drive velocity-force curve by plotting the 
average steady state velocity (V) versus the magnitude 
of the pulse drive Fd- This is shown in Fig. 2(a) and 
Fig. 2(c) for driving along y and x, respectively. For 
driving in the easy or y-direction, the critical force F^ is 
lower than for driving in the hard or x-direction. 

There is a sharp jump to a higher value of (V y ) just 
above F% = 0.165 in Fig. 2(a). For F V D < 0.165, the 
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FIG. 2: (a) The average steady state velocity (V y ) vs F D for 
pulse driving in the y-direction. We use the same parameters 
as in Fig. 1 with F v = 0.9 and p p = 0.38. Here the drive is 
suddenly increased from zero to Fd and the system settles into 
a steady state after a transient time r. (b) The corresponding 
d{V y ) /dF D , with a sharp peak indicating the transition to a 
moving stripe phase for Fp > 0.165. (c) (V x ) vs Fg for 
the same system, (d) The corresponding d(V x }/dF^ vs Fg 
showing a sharp peak at the transition to a perpendicularly 
translating stripe state. 



pinned stripe state undergoes plastic distortion when it 
moves, but for F^ > 0.165, the pinned stripe is able 
to depin directly into an ordered moving stripe state, 
resulting in the jump in mobility. An example of the 
ordered motion is shown in Fig. 3(e) and Fig. 3(f) for 
= 0.2. For Fjj < 0.165, the stripe structures are par- 
tially destroyed when plastic flow occurs at depinning. 
For weakly plastic flow, the moving stripes persist tran- 
siently for a period of time before breaking apart and 
rcpinning, as in Fig. 3(a) and Fig. 3(b) at F^ = 0.05. In 
contrast, the stripes are disordered in the strongly fluc- 
tuating plastic flow phase when a portion of the particles 
are pinned while other particles are mobile, as illustrated 
in Fig. 3(c) and Fig. 3(d) for F D = 0.125. The transition 
from the disordered plastic flow regime to the ordered 
moving stripe regime appears as a pronounced peak in 
d{V y )/dF v D at F V D = 0.165, as shown in Fig. 2(b). 

For the strongly fluctuating plastic flow regime found 



for 0.1 < Fp < 0.165, the initial transient motion con- 
sists of an elastically moving stripe state in which all the 
particles are moving. To characterize the time r required 
before the system reaches a steady state after application 
of a pulse drive, we analyze the time series of V y such as 
those plotted in Fig. 4 for F V D = 0.11, 0.125, 0.15, and 
0.17, where each point is averaged over 200 simulation 
time steps. The system starts with a higher value of V y 
which persists for a time that increases with increasing 
Fp before dropping to the lower steady state value of V y . 
The initial motion associated with the higher value of V y 
is a metastable moving ordered stripe. After the tran- 
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FIG. 3: The particle positions (black dots) and particle 
trajectories (black lines) for the system in Fig. 2(a,b) for 
pulse driving in the y-direction. The particle trajectories 
are traced over equal times in panels (a), (c), and (e). (a) 
At F^ — 0.05 there is a filamentary flow pattern along the 
stripes, (b) A snapshot of only the particle positions from (a) 
shows that the stripe structure is partially preserved, (c) At 
Fp = 0.125, fluctuating plastic flow occurs in which channels 
of moving particles intertwine and mix while other particles 
remain pinned, (d) The particle positions only for the system 
in (c) indicate that the stripe structures are completely disor- 
dered, (e) At F D = 0.2, above the peak in d{V y ) / dF D shown 
in Fig. 2(b), all the particles are moving in an ordered stripe 
phase, (f) A plot of only the particle positions from (e) shows 
the ordered stripes. 



sient time t, the stripe breaks apart and a portion of the 
particles become pinned, producing the drop to the lower 
value of V y . The particle positions shown in Fig. 3(d) at 
F^ = 0.125 are illustrated at a point in time after the 
transient ordered stripe state broke apart. 

For F D > 0.165, the system remains in the ordered 
moving stripe state within the entire simulation time win- 
dow, which includes simulations ten times longer than 
shown in Fig. 4. For F^ < 0.165, the transient time r 
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FIG. 4: The velocity V y in the indirection averaged over 200 
simulation time steps vs simulation time for pulse drives of 
F% = 0.11, 0.125, 0.15, and 0.17, from bottom to top. In 
the lower three curves, which fall in the strongly fluctuating 
plastic flow regime, the transition from a higher to a lower ve- 
locity occurs when the ordered stripe structure breaks apart, 
such as when the structure seen in Fig. 3(f) turns into a frag- 
mented structure of the type shown in Fig. 3(d). The upper 
curve is at a drive above the transition to the moving ordered 
stripe regime. 



during which the metastable ordered moving stripe ex- 
ists increases with increasing F^ . It is possible that after 
extremely long times, even for Fd > 0.165 the stripe 
state could break apart, resulting in a shift of the peak 
in d(Vy)/dF^ to higher Fp. Figure 2(b) shows that there 
is a linear increase in d(V y ) jdF^ for 0.1 < Fp < 0.16, be- 
low the large peak. Within this range of Fp , the steady- 
state flow is strongly fluctuating as shown in Fig. 3(c) 
and Fig. 3(d). For 0.04 < F% < 0.1, the (V y ) versus F D 
curve increases very slowly above the depinning transi- 
tion, as also indicated by the small value of dF^/dFo 
in Fig. 2(b). In this range of Fg\ the flow is still plas- 
tic and a portion of the particles remain immobile while 
others flow past; however, the character of the plastic 
flow differs from the strongly fluctuating flow found for 
0.1 < Fjj < 0.16. The low drive plastic flow takes the 
form of filamentary flow along the stripes, as illustrated 
in Fig. 3(a) for F^ — 0.05 where one river of particles 
flows along the stripe. Figure 3(b) shows that the parti- 
cles still retain much of the stripe structure, in contrast 
to the strongly fluctuating flow illustrated in Fig. 3(d) 
where the stripe structure is nearly lost. Another differ- 
ence is that the plastic flow in Fig. 3(c) involves a signif- 
icant transverse diffusion of particles in the x-direction, 
so that over time the particles can mix throughout the 
system. For the low drive plastic flow in Fig. 3(a), there 
are some early time particle jumps transverse to the drive 
from one stripe to another; however, these events vanish 
in the long time limit and there is no steady state dif- 
fusion in the cc-direction, even though the shape of the 
filamentary flow within one stripe may change slightly 
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FIG. 5: A combined plot of (V x ) vs Fg (dark line) and (V y ) 
vs (light line) from Fig. 2(a) and Fig. 2(c) highlights the 
difference in the depinning thresholds F§ and Fq and the 
transport characteristics. 



over time. 

In Fig. 2(c) we plot (14) versus Fg for the same system 
in Fig. 2(a), while the corresponding d(V x )/dF'f ) appears 
in Fig. 2(d). For the x-direction pulse drive, we do not 
find any regime where the stripes can dynamically re- 
orient and align themselves in the x direction. Instead, 
the system passes directly into a moving stripe phase in 
which the stripe orientation remains perpendicular to the 
direction of stripe motion. The transition into the sliding 
stripe phase occurs at the peak in d(y x )/dF'f ) shown in 
Fig.2(d) at Fg = 0.224. This peak falls at a higher value 
of Fg than the value of Fp of the peak in Fig. 2(b). This 
is because stripes moving perpendicular to their orienta- 
tion are much more susceptible to breaking apart than 
stripes moving parallel to their orientation. A higher 
pulse driving force reduces the effectiveness of the pin- 
ning and reduces the tendency for plastic flow, so a stripe 
moving perpendicular to its orientation is stabilized at a 
higher drive than a stripe moving parallel to its orienta- 
tion for the same pinning strength. 

Figure 2(d) also shows a smaller maximum in 
d{V x )/dFf ) near Fg = 0.135, corresponding to a change 
in the character of the plastic flow. For 0.135 < Fg < 
0.224, the steady state flow is dominated by strong plas- 
tic rearrangements where the stripe structure completely 
breaks apart. In contrast, for 0.06 < F'fy < 0.135, the 
flow is more filamentary and is composed of a small num- 
ber of slowly changing channels. For 0.17 < Fg < 0.22, 
we also find a metastable effect similar to that shown in 
Fig. 4 for driving in the y-direction. An ordered stripe 
can slide perpendicular to the direction of the stripe ori- 
entation for a period of time which increases with in- 
creasing Fg before the stripe breaks apart. 

In Fig. 5 we plot (14) versus Fg on the same panel 
as (V y ) versus Fq to highlight the regions of anisotropic 
transport. For 0.15 < Fg' y < 0.224, the velocity is higher 



for driving along the easy y direction. In contrast, for 
0.08 < Fp y < 0.15, (V y ) falls below (V x ). This occurs be- 
cause the pulse drive measurements preserve some of the 
initial structure of the oriented stripes. In the low drive 
regime 0.08 < F^' v < 0.15, much of the flow is filamen- 
tary. For driving along the stripe, Fp, the filamentary 
flow settles quickly into a few nonfluctuating channels, 
while for driving against the stripe, Fp, the filamentary 
flow forms fluctuating plastic channels which generally 
have large velocity pulses, producing a larger average ve- 
locity (V x ). For 0.025 < F^' y < 0.04, anisotropy appears 
due to the differing critical depinning forces in the two 
directions; here, flow only occurs in the y-direction but 
is absent in the x-direction. 



A. Transient Times 

We next analyze in detail the time required for the sys- 
tem to achieve steady state flow under a pulse drive. This 
time grows rapidly near the transition between the disor- 
dered and ordered flow states, as shown in Fig. 4. Recent 
experiments and simulations of periodically sheared par- 
ticle assemblies under suddenly applied shear have shown 
evidence for a diverging time to reach a steady state upon 
approaching a dynamic phase transition 5 - 3 .. Recent simu- 
lations of the plastic depinning of repulsively interacting 
particles also revealed that the transient time to reach 
a steady state under a pulsed drive diverges as a power 
law as the depinning transition is approached from either 
side 5 ^. This suggests that analyzing the transient time 
required to reach a steady state can be a useful diagnostic 
for probing changes in dynamical states. 

In Fig. 6(a) we plot the transient time r taken by the 
system in Fig. 2(a) to reach a steady state after the appli- 
cation of a pulse drive in the y-direction. For Fjj < 0.04, 
below the critical depinning force Fg, there is still a finite 
transient time during which the system organizes into a 
pinned state. The time required to reach the pinned state 
increases as F% is approached from below, while the time 
required to reach a steady moving state increases as F% is 
approached from above. For Fjj > 0.165, when the sys- 
tem passes directly into a moving ordered stripe state, 
the steady state velocity is reached very quickly and r is 
very small. In the region 0.1 < Fp < 0.165, the system 
depins into a metastable moving ordered stripe state that 
breaks apart after a time r as shown in Fig. 4. As Fp 
approaches F% = 0.165 from below, r increases since it 
takes increasingly longer times to trigger the instability 
that results in the fragmentation of the stripes. The rapid 
increase of r suggests that r may diverge at the transi- 
tion to the moving stripe phase; however, our results are 
not accurate enough to establish whether this divergence 
has a power law form. We find that the peak in t at 
Fp = 0.165 is asymmetric. In comparison, the dynamic 
phase transitions studied in the shearing systems pro- 
duced symmetric diverging time scales on both sides of 
the transitio n 53 ' 54 . Fig. 6(a) also shows a peak in r cen- 
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FIG. 6: The time r required for the system to reach a steady 
state velocity after an applied pulse drive for the system in 
Fig. 2. (a) r vs for driving in the y-direction. The first 
peak at — 0.04 corresponds to the depinning transition, 
the second peak centered at F^ = 0.1 appears at the tran- 
sition from the filamentary flow illustrated in Fig. 3(a,b) to 
the fluctuating flow state shown in Fig. 3(c,d), and the third 
peak at F^ — 0.165 is associated with the transition to a 
moving ordered stripe state, (b) r vs Ff> for driving in the 
x-direction. 



tered near F^ = 0.09, which corresponds to the location 
of the change in slope of (Vy) versus F^ in Fig. 2(a). At 
this drive, there is a change from the filamentary plastic 
flow channels shown in Fig. 3(a,b) to the rapidly fluctuat- 
ing disordered plastic flow channels shown in Fig. 3(c,d). 
The fact that r also increases in this region is further ev- 
idence that there can be dynamical phase changes even 
within the plastic flow regime. Finally, there is another 
peak in r near F^ — 0.04 at the depinning transition. 
These results show that peaks in the transient time can 
be used to detect changes in the flow characteristics of 
these systems. 

In Fig. 6(b) we plot t for pulse driving in the in- 
direction. Near Fg = 0.225 there is a peak in t as- 
sociated with the transition to the moving perpendicu- 
lar stripe phase. Within the range 0.21 < Ff, < 0.225, 
where the value of r is locally enhanced, the system forms 
a metastable state of stripes moving perpendicular to 
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time x 200 



FIG. 7: Solid line: V y vs time for Fg = 0.085 for the sys- 
tem from Fig. 5, showing the decay of the system into the 
filamentary plastic flow state. Dashed line: a power law fit to 
V v {t) oc t~ a with a = 0.35. 



their length. Eventually, the stripes break apart, and the 
time required for this to occur increases with increasing 
Ffy until for high enough Ffj the moving perpendicular 
stripe structure becomes stable rather than metastablc 
and r drops back to a small value. The peak in r at 
Ffj = 0.225 in Fig. 6(b) is similar to the peak found for 
driving in the y-direction in Fig. 6(a) near F^ = 0.165. 
Figure 6(b) shows that there is another peak in r centered 
at Ff) = 0.06 corresponding to the transition from fila- 
mentary plastic flow to strongly disordered plastic flow. 
Overall, the transient times for driving in the x-direction 
are smaller than for driving in the y-direction. For driv- 
ing along the y-direction, plastic flow channels can form 
which do not distort the aligned stripe pattern, permit- 
ting the system to remain in a metastable state for longer 
times before falling into the disordered steady state stripe 
structure. For driving along the x-direction, the stripe 
structure is more strongly disordered even at lower drives, 
so the system is closer to the disordered steady state 
stripe structure from the beginning and spends a shorter 
amount of time in the metastable state. 

The filamentary plastic flow regime that appears for 
0.04 < Fjj < 0.1 is associated with large values of r as 
shown in Fig. 6(a). In this regime, the time decay of 
V y to its steady state value differs from the decay in the 
strongly fluctuating plastic flow regime shown in Fig. 4, 
where V y remained roughly constant before dropping rel- 
atively rapidly to a lower value. Instead, in the filamen- 
tary regime V y follows a continuous stretched exponential 
form or power law, as illustrated in Fig. 7 for F^ — 0.085. 
The dashed line is a power law fit to V y oc t~ a with 
a = 0.35. We find an equally good fit of the curve to a 
stretched exponential form. The exponents from the fits 
do not appear to be universal and change when we take 
measurements from the other plastic flow regimes. Our 
results indicate that within the filamentary plastic flow 



phases, very long transient times can occur. 



The order-disorder transition at F^ — 0.165 between 
the lower drive strongly fluctuating plastic flow phase 
shown in Fig. 3(d) and the higher drive moving ordered 
stripe phase shown in Fig. 3(f) exhibits metastability and 
has a diverging time scale only on the low drive side of 
the transition, as indicated in Fig. 6(a). These features 
strongly suggest that this transition is first order in na- 
ture and that the details of the transition are strongly 
affected by the initial conditions of the moving state. For 
example, it is possible to obtain a reversed metastabil- 
ity by starting the system in a disordered configuration 
and applying a pulse drive Ff-> > 0.165. In this case, the 
moving system remains disordered and travels at a lower 
velocity until an instability causes the stripe structure to 
form with a corresponding increase in the velocity. This 
reversed metastability shows diverging transient times as 
the transition is approached from above, but has no di- 
verging time scales when the transition is approached 
from below. The metastability of the ordered and dis- 
ordered states resembles the superheating or supercool- 
ing recently observed for systems with first order phase 
transitions. Very similar dynamical superheating and 
supercooling effects were found in superconducting vor- 
tex systems in experiments^ and three-dimensional (3D) 
simulations^.. The vortex system undergoes a disorder- 
induced first order phase transition, and the effective dis- 
order changes when the system is prepared in different 
states. Computational studies of 2D vortex systems in- 
teracting with disorder have shown that there is either a 
continuous order to disorder transition or a crossover, so 
hysteresis, superheating and supercooling do not appear. 
For many stripe forming systems in two dimensions, tran- 
sitions from ordered to disordered states in equilibrium 
and in the absence of quenched disorder are first order 
in nature^. Our results suggest that the first order na- 
ture of the equilibrium transitions persists for some of 
the transitions between nonequilibrium states. 



At the transition between the low drive filamentary 
plastic flow phase and the higher drive strongly fluctu- 
ating plastic flow phase, Fig. 6(a) indicates that there 
are diverging transient times on both sides of the transi- 
tion. This behavior is similar to the diverging transient 
times found for 2D plastic depinning. There is evidence 
that the plastic depinning is an absorbing phase tran- 
sition falling in the directed percolation class^. Deter- 
mining whether the depinning of the stripe system or 
the filamentary plastic flow to strongly fluctuating plas- 
tic flow transition are also nonequilibrium phase transi- 
tions falling in the directed percolation class is beyond 
the scope of this work; however, our results suggest that 
the stripe system may be an ideal system in which to 
examine the nature of nonequilibrium transitions since it 
exhibits several different types of flow phases. 
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FIG. 8: (a) (V y ) vs for a pulse drive system with p p = 0.38 
and a lower pinning force of F p = 0.125. A two step depin- 
ning occurs, with uncoupled stripes depinning initially fol- 
lowed by the depinning of coupled stripes, (b) The corre- 
sponding d(V y )/dFp has a double peak indicating the two 
step depinning process, (c) (V x ) vs Fp for the same system, 
(d) The corresponding d{V x )/dFp shows a single step elastic 
depinning of the stripes. 
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FIG. 9: A combined plot of (V x ) vs Ff, (dark line) and (V y ) vs 
F% (light line) for the system in Fig. 8 with F p = 0.125 high- 
lighting the transport anisotropy. Inset: The particle posi- 
tions (dots) and trajectories (lines) for the system in Fig. 8(a) 
in the decoupled stripe regime at F^ = 0.0025 where a por- 
tion of the stripes are moving while others are pinned. 



with j5 — 0.35. This is followed by a crossover to 
(V x ) oc at higher drives. The behavior agrees well 
with the depinning of a harmonic elastic string driven 
over a random substrate, where an exponent /3 = 0.33 
is observed^. For driving in the y-direction, as shown 
in Fig. 8(a), (V y ) versus Fd indicates that a two step 
depinning process occurs. The initial depinning involves 
the flow of individual stripes, while other stripes remain 
pinned or move at different velocities, as illustrated in 
the inset of Fig. 9 for Fp = 0.0025. At higher drives, 
the remaining stripes depin and become coupled to the 
other moving stripes, resulting in an elastic flow. The 
two peaks in d(V y )/dFp shown in Fig. 8(b) fall at the 
locations of the two depinning transitions. A similar 
type of two step, layered depinning transition was pre- 
dicted for anisotropic charge density wave (CDW) sys- 
tems, where CDWs first depin separately and flow inde- 
pendently from one layer to the next, and then recou- 
ple at higher drives^. Mean field models also predict 
that layered systems should show a coupling-decoupling 
transitio n 59 i 60 , while two-layer models predict that co- 
existence of moving and pinned phases should occur in 
2D systems^. In the main panel of Fig. 9, we plot 
(V x ) versus Fg and (V y ) versus F|, together in order to 
highlight the transport anisotropy which disappears for 
Fp V > 0.04 when fully elastic flow is established. The 
transient times r for the F p — 0.125 system are much 
shorter than those in the F p = 0.9 system, where plastic 
depinning occurred. We find an increase in r just below 
each depinning transition in the F p = 0.125 system for 
both x and y direction driving. For x-direction driving, 
there is a single peak in r below the depinning thresh- 
old which is associated with a small amount of rough- 
ening of the stripe structure that occurs just before de- 
pinning. Above each depinning transition in the weak 
pinning system, t is extremely small. For y-direction 
driving, the peak in r is broader within the sliding plas- 
tic flow phase found below the second depinning transi- 
tion. In systems with even weaker pinning, F p < 0.05, 
the depinning is elastic for both driving directions and 
the transport anisotropy is significantly reduced. 



IV. CONTINUOUS FORCE SWEEP 
MEASUREMENTS 

A. Strong Disorder 



B. Pulse Measurements for Weak Disorder 

We next consider pulse drive measurements for a sys- 
tem with the same pinning density p p = 0.38 but with a 
weaker disorder strength of F p = 0.125. In Fig. 8(c) we 
plot (V x ) versus Fg and in Fig. 8(d) we show the cor- 
responding d(V x )/dFp. For this value of F p , the stripe 
structure remains ordered. The single step depinning for 
x-direction driving is elastic and has (V x ) oc (Fg — F*)° , 



We now examine the case where the applied drive is 
slowly incrementally increased in a single sweep as op- 
posed to the sudden application of the drive discussed in 
the previous section. We use a driving force increment 
of AFo = 0.00025 applied every 25000 simulation time 
steps in a system with F p — 0.9 and p p — 0.38. In Fig. 10 
we plot the resulting (V x ) versus Fg and (V^) versus F^ 
together. Many of the transport features are the same 
as those shown for the pulse drive in Fig. 5, such as the 
lower depinning threshold for driving in the y-direction 
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FIG. 10: A combined plot of (V x ) vs F£ (dark line) and (V y ) 
vs (lighter line) for a system with F p = 0.9 and the same 
parameters as in Fig. 2 but with a continuously swept drive. 
There are several steps in (K,) that appear below the drive 
Fp V — 0.32 at which the two curves meet. Inset: A blow-up 
of the region near F^' y = 0.24 shows a crossing of the (V x ) 
and (V y ) curves, indicating that (V y ) has a smaller slope and 
that therefore fewer particles are moving for driving in the y- 
direction than for driving in the x-direction. The crossing of 
{V x ) and (Vy) does not occur for the pulse drive measurements 
shown in Fig. 5. 



y . y 



(a) x (b) x 

FIG. 11: (a) Image of particle positions in the reordered 
phase at Fp = 0.27 in the swept drive system from Fig. 10. 
In addition to the stripe structure, there are some scattered 
particles pinned between the stripes, (b) Image of particle 
positions in the reordered phase at Ff, = 0.27 in the same 
system. The stripes have reoriented in the x-direction and 
there are no pinned particles between the stripes. 

and the anisotropic flow centered near F^' v = 0.2. The 
transition to the flowing stripe state remains sharp for 
the continuous sweep drive in the y-direction, but for the 
x-direction driving the state in which the stripes flow 
perpendicular to their orientation is lost. Instead, the 
system reorders into a stripe state oriented along the 
^-direction, as shown in Fig. 11(b) for Fjj = 0.27. In 
the pulse drive system, the perpendicular moving stripe 
state can be stabilized by dynamical quenching, but for 
the swept drive, the system passes through an extensive 



plastic flow phase which destroys any memory of the ini- 
tial perpendicular stripe orientation. 

The inset of Fig. 10 highlights that under the swept 
drive, a crossing of (V y ) and (V x ) occurs near F^' v = 0.24. 
There is no such crossing for the pulse drive, as shown 
in Fig. 5. The crossing of the curves indicates that al- 
though the system is in the moving stripe state for both 
drive directions, the slope of (V y ) is smaller than that of 
(V x ) in this regime. This is because a fraction of the par- 
ticles remain pinned for the y-direction drive, while all of 
the particles are moving for the x-direction drive. When 
fewer particles are moving, the slope of the velocity-force 
curve is reduced. Particles are able to remain pinned 
for the y-direction drive because they can be captured in 
pinning sites that are sufficiently far away from the neigh- 
boring stripes that they experience only repulsion from 
the particles in the stripes, and are out of the range of the 
attractive part of the particle-particle interaction poten- 
tial. This does not happen for the cc-direction drive be- 
cause during the stripe reorientation process, all pinned 
particles are eventually swept up into a moving stripe. In 
Fig. 11(a), the reordered stripe phase for Fp — 0.27 con- 
tains pinned particles that sit between the stripes rather 
than flowing with the stripes. In contrast, in the mov- 
ing stripe phase for a y-direction pulse drive shown in 
Fig. 3(f), there are no pinned particles between the mov- 
ing stripes. In the a;-direction swept drive moving stripe 
phase, all of the particles are moving and there are no 
pinned particles between the moving stripes, as shown in 
Fig. 11(b) for Ffj = 0.27. The transition into the moving 
stripe phase for swept y-direction driving is rapid, as in- 
dicated by the jump in (V y ) at Fp = 0.185 in Fig. 10. For 
x-direction driving, the reordering transition is more con- 
tinuous, permitting more meandering of the stripe pat- 
tern during the stripe formation process. This allows all 
of the pinned particles to be attracted gradually into the 
moving stripe structure. 

Once the driving force becomes strong enough, the 
pinned particles surrounding the moving stripes for the 
y-direction swept drive depin and join the moving stripe 
structures. The depinning of the individual particles pro- 
duces the step features in (V y ) near Fp = 0.28 in Fig. 10. 
For drives above the depinning threshold of all of the 
pinned particles, the velocity response is isotropic, as 
shown in Fig. 10 for F^ v > 0.32. Our results indicate 
that even for pinning strengths strong enough to induce 
plastic flow and a subsequent reordering of the stripe 
structure, some memory of the initial ordering of the 
stripe phase is retained up to relatively large values of the 
driving force. We note that based on the extended tran- 
sient behavior found in the plastic flow regimes for the 
pulse drives, it is possible that if the swept drives were ap- 
plied with even slower drive increments, the anisotropic 
response could be lost at lower drives if the system is 
given more time to slowly mix in the plastic flow state. 
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FIG. 12: A combined plot of (V x ) vs (dark line) and 
(V y ) vs (lighter line) for a system with p v — 0.38, a lower 
F p = 0.225, and a swept drive. The depinning in the in- 
direction is accompanied by plastic distortions and the partial 
breaking apart of the stripe structure. Depinning in the y- 
direction occurs by the sliding of stripes past one another in 
a manner similar to that illustrated in the inset of Fig. 9. 
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B. Pinning Strength Dependence and Peak Effect 

When we perform swept drive measurements on a sys- 
tem with the weaker pinning strength of F p = 0.125, we 
obtain velocity-force curves that are nearly identical to 
those shown in Fig. 8 for a pulse drive measurement. We 
attribute this to the lack of plastic flow in the weakly 
pinned system. Without plastic flow, the stripe struc- 
ture never breaks apart and memory of the initial stripe 
orientation is never lost, so the same stripe orientation 
appears, regardless of whether the drive is slowly swept 
or suddenly applied. By increasing F p slightly, we reach 
a state where the x-direction depinning is plastic and ac- 
companied by the breaking apart of the stripe structure, 
while the y-direction depinning occurs by the sliding of 
some of the stripes, with the stripe structure maintained 
intact. This situation is illustrated in Fig. 12, where we 
plot (V x ) versus Ff) and (V y ) versus for a swept drive 
system with F p = 0.225. 

By conducting a series of simulations for varied F p , 
we map the anisotropy in the depinning thresholds F* 
and i™, as shown in Fig. 13(a). For x-direction driv- 
ing, Fj? increases monotonically with increasing F p for 
< F p < 0.175. Within this range of F p , the depin- 
ning is elastic and the stripes move perpendicularly to 
their orientation. Each particle maintains its position 
in its original stripe. Just above F p — 0.175, there is a 
sudden increase in F* corresponding to the onset of plas- 
tic distortions of the stripe structure at depinning. For 
driving in the y-direction, F% continuously increases with 
increasing F p for < F p < 0.25. Over this range of F p , 
the stripes depin elastically and either all depin simul- 
taneously for F p < 0.05, or depin as individual sliding 



FIG. 13: (a) The critical depinning forces F" (filled circles) 
and F£ (open squares) for driving in the x and {/-directions, 
respectively, plotted vs F p . The sharp jumps in the depinning 
forces are associated with changes from elastic flow to plastic 
flow and from plastic flow to sliding ordered flow, (b) The 
corresponding dF" /dFp vs F p (dark line) and dF^ /dFp vs F v 
(light line) curves more clearly show the onset of the different 
depinning phases, (c) The anisotropy ratio R — F* / F^ vs F p 
shows that the change in anisotropy can be associated with 
different depinning regimes. 



stripes for 0.05 < F p < 0.25 via the mechanism illus- 
trated in the inset of Fig. 9. For F p < 0.05, where our 
system behaves elastically, we find that F-f increases ap- 
proximately quadratically with F p which is the expected 
behavior for elastic depinning22,. Just above F p = 0.25 
there is a sharp increase in Ffj when the stripes begin to 
depin plastically along individual stripes in the manner 
shown in Fig. 3(a). 

For F p > 0.4, both i™ and F* begin to saturate, as 
shown by the anisotropy ratio R = F^/Fy plotted in 
Fig. 13(c). The saturation arises due to the fact that the 
number of pinning sites N p is only slightly higher than 
the number of particles N . As a result, when F p is large 
enough, the initial depinning is dominated by intersti- 
tially pinned particles. These particles are not trapped 
by one of the randomly located pins, but are instead held 
in place by interactions with neighboring pinned parti- 
cles. The critical force for the depinning of particles 
trapped by pins increases linearly with increasing F pi but 
the critical force for the depinning of interstitially pinned 
particles is determined only by the particle-particle inter- 
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action potential and is not altered by increasing F p . In a 
system with a much higher pinning density, every particle 
would be trapped by a pin and the depinning threshold 
would show the expected linear increase with increasing 
F p . 

The sharp increases in the depinning thresholds associ- 
ated with transitions from elastic to plastic flow or from 
weakly plastic flow to a more strongly plastic flow resem- 
ble the phenomenon observed for depinning of vortices 
in type-II superconductors, where a peak in the depin- 
ning threshold has been connected with the disordering 
of the vortex lattice^ 1 ' 55 ' 56 ! 63 ' 64 . In the disordered or 
plastically flowing systems, the particles can more read- 
ily adjust their positions to take advantage of the energy 
of a randomly located pinning site without paying the 
large energy cost required to distort an elastic or ordered 
particle lattice. In studies of 2D vortex systems, as the 
vortex lattice is softened the system becomes more disor- 
dered and the depinning threshold increases continuously, 
producing a peak effect that is continuous rather than 
sharp 65 . This may be related to the fact that 2D systems 
of particles with repulsive long range interactions lack 
first order melting or disordering transitions. Many 3D 
vortex systems have a first order transition from an or- 
dered vortex structure to a disordered one, and the peak 
effect phenomena observed in these systems is very sharp. 
Even through the stripe system described here is 2D, the 
pairwise particle interactions are not strictly repulsive. 
Previous studies in the absence of quenched disorder us- 
ing the same model produced results that suggest that 
the thermal melting of the stripe and clump systems is a 
first order transition. Other 2D studies of stripe-forming 
systems with competing interactions also find first order 
melting transitions for many of the phases^. These re- 
sults suggest that the first order peak effect phenomenon 
found for vortex systems may also generically occur for 
stripe and pattern forming systems in the presence of 
quenched disorder 

To further characterize the changes in F c we plot 
dF*/dF p and dFV/dFp vs F p in Fig. 13(b). There is 
a peak in dF*/dF p near F p = 0.185 corresponding to 
the transition from ordered stripe flow to partially plas- 
tic flow. A second, much broader peak appears near 
F p = 0.35 at the point where the stripes break apart 
completely. For 0.185 < F p < 0.35, the stripes driven 
in the ^-direction reorder into a perpendicularly moving 
stripe state at high Ffy in spite of the fact that weakly 
plastic flow occurs above depinning. For F p > 0.35, the 
plasticity at depinning becomes much stronger and the 
stripes reorder into a parallel moving state at high Ff). 
For driving in the y-direction we find a single sharp peak 
in dFy/dF p near F p = 0.25 corresponding to the appear- 
ance of plastic flow along the stripes. We note that the 
transition at F p — 0.075 into the pinned-sliding phase il- 
lustrated in the inset of Fig. 9 is not associated with any 
sharp features in dF^ /dF p . 

The appearance of different flow phases can also be de- 
tected in the plot of the anisotropy ratio R in Fig. 13(c). 
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FIG. 14: Dynamic phase diagrams showing the onset of de- 
pinning (lower heavy line) and the onset of the dynamic stripe 
reordering (upper light line), (a) F£, vs F p . The dashed line 
separates the moving perpendicular stripe phase from the dy- 
namically reordered parallel stripe phase. In order for the par- 
allel stripe phase to form, the system must first pass through 
a strongly fluctuating plastic flow phase, (b) vs F p . The 
dashed line indicates the change from sliding stripe plastic 
flow to plastic flow in which the stripe structure breaks apart. 



For example, the peak in R over the range 0.185 < F p < 
0.3 occurs when the flow for x-direction driving is plastic, 
while the flow for y-direction driving remains elastic in 
the stripe sliding state. The increase in R at F p = 0.06 
corresponds to the transition from elastic depinning in 
both directions to elastic depinning for x-direction driv- 
ing and individual stripe sliding for y-direction driving. 
For F p > 0.6, the anisotropy begins to saturate. 



C. Dynamic Phase Diagram 

By identifying features in the velocity force curves and 
associating them with different moving states, we con- 
struct a dynamic phase diagram for both x and y direc- 
tion driving for F p — 0.225. In Fig. 14(a) we plot the 
phase diagram for Fp vs F p , indicating the location of 
the depinning curve and the transition into an ordered 
stripe state. For F p < 0.185, the stripes depin elasti- 
cally into a perpendicular moving stripe state, while for 
0.185 < Fp < 0.35, the stripes depin with a small amount 
of plastic distortion into the same perpendicular moving 
stripe state. The line in Fig. 14(a) marking the transition 
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from plastic flow to reoriented ordered parallel moving 
stripes increases roughly linearly with F p for F p > 0.4. 

In Fig. 14(b) we plot the phase diagram of versus 
F p for driving in the y-direction. In this case the stripes 
are always oriented in the direction of the drive. The 
range of the plastic flow regime grows with increasing F p 
and in general the onset of the moving stripe phase occurs 
at lower drives than those at which the parallel moving 
stripes form for x-direction driving. The small dashed 
line indicates the transition from the plastic flow in which 
the stripe structure is destroyed for F p > 0.275 to the 
state where moving stripes slide past pinned stripes for 
0.075 < F p < 0.275. 

The dynamic phase diagram for the stripe system con- 
tains a larger number of phases than dynamic phase dia- 
grams observed in systems with purely repulsive particle- 
particle interactions moving over random disorder. For 
example, in 2D vortex systems the dynamic phases con- 
sist only of a pinned state, a plastic flow state, and a 
moving partially ordered state, and the transitions be- 
tween these states are continuous. In the partially or- 
dered moving state, the particles are not fully crystal- 
lized but develop a smectic type of ordering and flow in 
evenly spaced channels aligned with the direction of the 
drive. The channels of flow may be coupled or partially 
coupled^r— . This is similar to the state we observe in 
which the stripes reorient in the direction of drive for 
sufficiently high drive and sufficiently strong pinning. In 
the stripe system, the stripe reordering transition is more 
consistent with a first order phase transition rather than 
the continuous or crossover behavior found for particles 
with purely repulsive interactions. Studies of driven sys- 
tems with quenched disorder, where the particle-particle 
interactions are more complicated than the purely repul- 
sive case, have shown that it is possible to have a coexis- 
tence of different moving phases, which is consistent with 
having first order phase transitions between the moving 
phases^. 




FIG. 15: (a) F£ (filled circles) and F? (open squares) vs 
p p for a system with F p — 0.225 and R p — 0.3. (b) The 
corresponding R = F% /F% vs p p shows that at large p p the 
anisotropy is reduced, (c) F% (filled circles) and F* (open 
squares) vs p p for a system with F p = 0.7 and R p = 0.3. (d) 
The corresponding Rvsp p . Here the anisotropy drops nearly 
to 7? = 1 at high p p 
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FIG. 16: (a) F* (filled circles) and F c y (open squares) vs R p 
for a system with F p — 0.225 and p p — 0.38. (b) The ratio 
R = F® /F* vs R p showing that the anisotropy increases with 
increasing R p and saturates at high R p . 



D. Changing Disorder Density and Radius 

In Fig. 15(a) we plot F* and F-f versus p p for a sys- 
tem with F p = 0.225 and R p = 0.3, and in Fig. 15(b) we 
show the resulting anisotropy ratio R — F* j ' F^ versus 
Pp. Here Fj? and F^ both increase monotonically with 
increasing p p . At low p p the anisotropy R is strongly re- 
duced and the two depinning curves come together when 
the depinning becomes elastic for both x and y-direction 
drives. A similar effect appears in Fig. 13(c) for low F p , 
where the onset of elastic depinning for both driving di- 
rections results in a reduced value of R. The anisotropy 
in Fig. 15(b) passes through a maximum near p p = 0.2 
before gradually falling back to 1.0 for increasing p p . The 
pinning density plays a more important role in deter- 
mining the anisotropy of the depinning for pre-formed 
stripe states than the pinning strength. When the av- 
erage distance between pinning sites is greater than the 



inter-stripe distance, depinning for y-direction driving in 
a system with strong pinning occurs via a combination 
of plastic flow of particles along some stripes while other 
stripes remain completely pinned. As a result, i™ is gen- 
erally lower than F* for the low density, strong pinning 
limit. This is shown in Fig. 15(c) where we plot i™ and 
F* versus p p for F p — 0.7. For higher values of p p , the 
plastic depinning along the stripes for y-direction driving 
is suppressed and the behavior becomes more isotropic. 
For systems with small p p and F p > 0.35, the anisotropy 
persists down to much lower values of p p than shown in 
Fig. 15(a); however, for p p > 1.0 the anisotropy vanishes 
completely. 

We next consider a system with F p = 0.225 and 
p p = 0.38 with varied pinning radius R p . We plot F* 
and F£ versus R p in Fig. 16(a), and the anisotropy ratio 
R = F*jpy in Fig. 16(b). Both F* and i™ increase 
with increasing R p . For low R p , F* increases faster 
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FIG. 17: (a) The time series of V y for a system with F p = 0.9, 
p p — 0.36, R p = 0.2, and = 0.04. For these parameters, 
the particles move in stationary plastic filaments, producing 
a periodic time-of-flight signal, (b) The power spectrum S(f) 
of the time series in (a) showing the characteristic peaks from 
the periodic signal. Here, / is in units of inverse simulation 
time steps, (c) Time series of V y for the same system at 
= 0.12 in the strongly fluctuating plastic flow regime, (d) 
The corresponding S(f) has a 1// feature at low frequencies 
as indicated by the dashed line. 

with increasing R p than i™ does, and for high i? p , the 
anisotropy R saturates. 

V. VELOCITY FLUCTUATIONS 

The dynamic phases can also be characterized by mea- 
suring the velocity fluctuations. Experiments on stripe 
forming systems previously demonstrated that narrow 
band noise, characterized by a periodic noise signal, and 
broad band noise, which lacks any characteristic frequen- 
cies, can occur and showed evidence for transitions be- 
tween the different types of nois o 31 ' 67 . Previous simu- 
lations of clump and stripe forming systems showed the 
presence of a 1// noise characteristic in the nonlinear 
portion of the velocity-force curve associated with the 
fluctuating plastic flow phase. At high drives where the 
stripes or clumps reorder, the noise becomes white with 
a weak narrow band or washboard frequency similar to 
that observed for a driven vortex lattice in the dynami- 
cally reordered regime^. 

Here we show that the stripe system exhibits many 
additional noise features near depinning. For y-direction 
driving near depinning, there can be filamentary plas- 
tic flow channels along the stripe with no particle diffu- 
sion from stripe to stripe, as shown in Fig. 3(a). Within 
this filamentary plastic flow regime, it is possible for the 
particle flow to be limited to one or a small number of 
individual winding channels which do not change over 
time and which have a characteristic time-of-flight for 
crossing the sample. This results in a periodic velocity 



signal such as that shown in Fig. 17(a) for a system with 
F p = 0.9 and = 0.04. Similar periodic filamentary 
plastic motion has been observed in vortex simulations 
performed just at depinning. Here one or two stable 
channels of moving particles form while the rest of the 
particles are immobile 68 . Evidence for filamentary flow 
has also been found in vortex experiments, where a se- 
ries of jumps and dips in the current- voltage curve were 
interpreted as indicating the opening of individual chan- 
nels of vortex flowfi. There are also several simulations of 
vortex systems showing transitions from narrow band fil- 
amentary flow to chaotic flow as the drive is increased 70 . 
The power spectrum S(f) of the velocity time series in 
Fig. 17(a) is shown in Fig. 17(b), and has characteristic 
narrow band noise peaks produced by the time-of-flight 
signature. As we increase F^ and permit the system to 
settle into a steady state, we observe a seres of transi- 
tions from ordered flows to fluctuating flows with broad 
band noise signatures. This behavior very similar to that 
found in 2DEG transport measurements 3 -^. We observe 
filamentary plastic flow from the depinning transition up 
to F V D = 0.1. Above this drive, the system transitions 
into the strongly fluctuating plastic flow regime shown in 
Fig. 3(c) in which the stripe structure is destroyed. In 
Fig. 17(c) we plot the time series of V y for the same sys- 
tem in Fig. 17(a) at a drive of = 0.12 in the strongly 
fluctuating plastic flow regime. The corresponding S(f) 
appears in Fig. 17(d). Here, V y (t) is strongly fluctuat- 
ing since the number of pinned particles in the system is 
continuously changing, and the power spectrum shows a 
1// noise characteristic at low frequencies. As the drive 
is further increased, the system reorders into a moving 
stripe state and the low frequency 1// spectral signal is 
lost. It is replaced by a weak narrow band noise signal 
similar to that observed in previous simulations for stripe 
reordering^. 

For driving in the a;-direction at F p = 0.9, the filamen- 
tary motion along the stripes of the type seen for driving 
along the y-direction is strongly suppressed. Instead, the 
strongly fluctuating plastic flow states with 1// noise sig- 
natures are more prevalent. When the stripes reorder by 
reorienting into the rr-direction, we find the same weak 
narrow band noise feature observed above the stripe re- 
ordering transition for driving in the y-direction. 

Near depinning in a system with weaker pinning of 
F p = 0.125, the sliding stripe phase produces a periodic 
signal as shown in Fig. 18(a) for i 7 ^ = 0.0025. The peri- 
odic signal also appears in the corresponding power spec- 
trum S(f) shown in Fig. 18(b). In this regime, the stripes 
are decoupled, so individual stripes are moving at slightly 
different velocities. This results in a more complex ve- 
locity signal composed of several similar frequencies. For 
higher drives, the moving stripes couple and the noise is 
more characteristic of a single periodic signal. In larger 
systems, near depinning there could be larger numbers 
of frequencies present since there are a larger number of 
stripes which can each move at different velocities. This 
would broaden the power spectrum; however, at larger 
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time 



FIG. 18: (a) The time series V y for a system with F p — 0.125 
and — 0.0025, in the sliding phase, (b) The corresponding 
power spectrum S(f). A strong narrow band noise signal 
appears, (c) Time series V x for the same system under driving 
in the ^-direction at Fg = 0.0375, just above depinning. (d) 
The corresponding S(f). The velocity signal is not periodic; 
however, the power spectrum in (d) has a Lorentzian shape 
with a flat spectrum at lower frequencies and a l// 2 spectra 
at higher frequencies, indicated by the dashed line. For high 
values of Ffj, the power spectrum develops a narrow band 
noise characteristic. 

drives, where the stripes couple, a strong narrow band 
noise signature should appear. 

In the weak pinning system of F p = 0.125, ^-direction 
driving produces elastic depinning. Just above the de- 
pinning transition, the noise signal is not periodic as 
shown by the plot of V x in Fig. 18(c) for Fg = 0.0375. 
There is, however, a characteristic noise frequency, as 
shown by the Lorentzian shape of S(f) in Fig. 18(d). At 
higher drives, the spectrum broadens and the Lorentzian 
peak frequency shifts to higher frequency with increasing 
drive. 



VI. THERMAL EFFECTS 

We next consider the effects of thermal fluctuations. 
In previous work for a stripe-forming system at the same 
density considered here but with no quenched disorder, 
we used diffusion and specific heat measurements to iden- 
tify a well-defined disordering temperature T m above 
which the stripe structures were completely destroyed^. 
Below T m , there was liquid- like particle motion along the 
length of the stripes, but the system behaved like a solid 
in the direction perpendicular to the stripes. This sug- 
gests that in the presence of quenched disorder, the stripe 
system might show considerable creep in the easy flow di- 
rection, but not in the hard direction for a finite tempera- 
ture at applied drives below the zero temperature depin- 
ning thresholds. To examine this, we consider a system 
with p p = 0.36 and F p = 0.7. At drives of Fg = 0.09 




FIG. 19: (a) {V x )/V (filled circles) and (V y )/Vo (open 
squares) vs T/T m for a system with p p — 0.36, R p = 0.2, 
and F p = 0.7 at applied drives of Ff, = 0.09 and F y D = 0.09, 
respectively. Here T m is the temperature at which the stripes 
melt in the absence of pinning and Vo is the velocity at 
which the particles would move in the absence of pinning. 
Anisotropic transport occurs for T/T m < 1.0. (b) Velocity 
anisotropy ratio {Vy)/{V X ) vs T/T m for the same system shows 
a diverging anisotropy at low temperatures. 



or Fjj = 0.09, the system is pinned in both directions at 
T = 0.0. In Fig. 19(a) we plot (V y )/V and (V x ) /V ver- 
sus T/T m . Here Vo is the velocity at which the particles 
would move in the absence of pinning. For T/T m < 0.25, 
there is almost no creep for x-direction driving but there 
is considerable creep for y-direction driving. As a result, 
at low temperatures the anisotropy of the velocity re- 
sponse diverges, as indicated by the plot of (V y )/{V X ) in 
Fig. 19(b). 

In Fig. 20(a) we plot the particle positions and trajec- 
tories for this system with Fg = 0.09 at T/T m = 0.25. 
Here there is liquid-like motion of the particles along 
the stripes but there is no diffusion perpendicular to the 
stripes. The stripe structure remains ordered although 
plastic creep is occurring. In the creep process, some 
particles along the stripe remain pinned while other par- 
ticles move around them along the length of the stripe. 
For 0.25 < T/T m < 1.0, creep occurs for both directions 
of drive but there is a larger amount of creep for driving 
in the y-direction. The anisotropy of the creep gradually 
diminishes as T approaches T m . In Fig. 20(b) we plot 
the particle positions and trajectories for = 0.09 and 
T/T m — 0.76. The stripe structure is still present but 
some hopping of particles from stripe to stripe occurs in 
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FIG. 20: The particle positions (dots) and trajectories (lines) 
at different temperatures for the system in Fig. 19 with Fjj = 
0.09. (a) At T/T m = 0.25, there is creep only along the 
{/-direction and no creep in the a;-direction. (b) At T/T m = 
0.76, the stripe structure is still present. There is considerable 
diffusion along the stripes and a smaller amount of diffusion 
across the stripes, (c) At T/T m — 1.266, the motion is no 
longer anisotropic, as indicated in Fig. 19, and the diffusion 
is isotropic, (d) The particle positions only from panel (c) 
show that the stripe structure is destroyed. 

the x-direction. For T/T m > 1.0, the creep anisotropy 
vanishes. Here the stripe structure is completely disor- 
dered and there is diffusion through the entire sample 
as shown in Fig. 20(c) and Fig. 20(d) for T/T m = 1.266. 
There is still some temporary trapping of particles by the 
pinning sites, indicated by the fact that (V Xt y)/Vo < 1.0. 
For higher temperatures, (V Xty )/Vo gradually approaches 
1.0 as the effectiveness of the pinning is diminished. In 
general we find that the creep anisotropy persists longer 
at lower pinning densities and that at higher pinning den- 
sities the creep anisotropy disappears. 

A. Thermally Induced Ordering 

We find an interesting effect in which thermal noise in- 
duces the formation of stripe order. At low temperature 
and for sufficiently strong disorder, the stripe structures 
are fragmented and destroyed. As the temperature is in- 
creased, it is possible for the thermal fluctuations to wash 
out the effectiveness of the pinning before the melting 
temperature of the stripe structure is reached. The re- 
sult is a floating ordered stripe. A similar effect has been 
observed for two-dimensional vorte»2i and colloid^ sys- 




FIG. 21: Phase diagram of T/T m vs F p for a system with 
p p = 0.38, R p — 0.2, and a {/-direction drive of = 0.22. 
Here T m is the melting temperature for the system without 
quenched disorder. The lower solid line separates the plastic 
flow phase from the moving stripe phase and the upper solid 
line separates the moving stripe phase from the moving liq- 
uid phase. The dashed line separates the plastic flow phase 
from the moving liquid phase. The phase diagram shows that 
for intermediate pinning strength, increasing the temperature 
can produce a transition from a disordered plastic flow phase 
into an ordered moving stripe phase. 
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FIG. 22: The particle positions (dots) for the system in 
Fig. 21 at F p — 0.9 for different temperatures. (a) At 
T/Tm = 0.1, the system is undergoing plastic flow and the 
particles are disordered, (b) At T/T m = 0.75, the tempera- 
ture has reduced the effectiveness of the pinning, permitting 
the formation of a moving stripe phase, (c) At T/T m = 1.3, 
the temperature is large enough to melt the stripe structure. 



terns interacting with periodic and random substrates. 
In the floating solid transition found in these studies, the 
vortices or colloids are pinned to the substrate at low 
temperatures, while at higher temperatures they float 
free of the substrate and form a triangular lattice. At 
still higher temperatures, the lattice disorders thermally. 
When the substrate pinning is strong enough, the float- 
ing solid phase disappears and the system passes directly 
from a pinned solid to a liquid state. 

To illustrate the formation of a floating stripe phase in 
our system in the presence of random disorder, in Fig. 21 
we plot a phase diagram of T/T m versus F p for a system 
with p p = 0.38, R p = 0.2, and F% = 0.22. Here T m is 
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the melting temperature of the stripes in the absence of 
quenched disorder. For low temperature and weak disor- 
der of F p < 0.3, the system is in a moving stripe phase. 
For 0.3 < F p < 1.75, at low temperatures the system is 
in the strongly disordered plastic flow state illustrated in 
Fig. 22(a) for F p = 0.9 and T/T m = 0.1. As the tempera- 
ture is increased, the effectiveness of the quenched disor- 
der is thermally destroyed and the system organizes into 
a moving stripe state aligned with the direction of drive, 
such as that shown in Fig. 22(b) for T/T m = 0.75. The 
system melts into a liquid for T/T m > 1.0 as shown in 
Fig. 22(c) for T/T m = 1.3. The phase diagram in Fig. 21 
indicates where these three phases occur. For F p > 1.6 
the system goes directly from a partially pinned plastic 
flow phase to an unpinned disordered liquid phase, as 
indicated by the dashed line. This result suggests that 
in systems with strong disorder, it is possible that the 
anisotropy may be weak at low temperatures but could 
increase for intermediate temperatures when the floating 
moving stripe structure forms. 

The ability of the system to form a floating stripe de- 
pends on the length scale of the quenched disorder as well 
as the disorder strength. If the disorder is composed of 
small well localized pins, as in our model, a floating solid 
phase is possible. In contrast, for long range quenched 
disorder the thermal noise will not be effective in wash- 
ing out the pinning and a floating solid phase will not 
occur. The presence of a floating solid phase makes it 
impossible to observe a thermally induced peak effect 
phenomena of the type found in superconducting vortex 
systems. In this peak effect, a thermally melted vor- 
tex lattice is softer and is able to better couple to the 
quenched disorder, increasing the depinning forced. If 
the thermal fluctuations destroy the effectiveness of the 
quenched disorder below the temperature at which the 
disorder-free equilibrium particle structure melts, a ther- 
mally induced peak effect can not occur. It is possible 
that adding long range correlations to the pinning would 
permit the appearance of a thermally induced peak effect 
in the stripe forming system. 



VII. SUMMARY 

We have examined the anisotropic dynamics of ori- 
ented stripes in a system with competing interactions. 
We focus on the regime where stripe structures form in 
equilibrium and in the absence of quenched disorder. Af- 
ter adding a random substrate, we drive the system par- 
allel or perpendicular to the original orientation of the 
stripes. We find anisotropic depinning thresholds and 
nonlinear velocity force curves. 

Under the sudden application of an external drive, the 
system settles into a steady state flow after a transient 
time that is determined by the structure of the steady 
state flow, such as a plastically flowing state or a moving 
ordered state. The transient time passes through peaks 
at the transitions between different dynamical states, 



such as from pinned to filamentary flow or from strongly 
fluctuating plastic flow to ordered flow. In addition to 
the stripe system, this type of transient measurement af- 
ter a sudden application of a drive could also be used 
in other driven systems such as superconducting vortices 
with quenched disorder, friction, and sliding charge den- 
sity waves. 

We observe different types of plastic flow which are 
determined by the direction of the drive relative to the 
stripe orientation. For driving parallel to the stripes, 
there is a phase in which the stripes remain ordered 
but are decoupled and can slide past one another. For 
stronger quenched disorder, plastic flow can occur within 
individual stripes while the overall stripe structure re- 
mains intact. In this case, the flow is filamentary and 
involves only a portion of the particles within the stripe. 
For stronger or denser quenched disorder, the stripes 
break apart and we find a strongly fluctuating plastic flow 
phase in which the transport properties are isotropic. For 
driving perpendicular to the stripe orientation, in addi- 
tion to plastic flow phases there can be elastic depinning 
of the stripes perpendicular to the drive for sufficiently 
weak disorder. 

As a function of disorder strength we find a sharp or- 
der to disorder transition in which the state above de- 
pinning changes from an ordered moving stripe structure 
to a plastic flow regime which tears apart the stripes. 
This order-disorder transition is accompanied by a sharp 
increase in the depinning threshold which is similar to 
the peak effect phenomenon observed near order-disorder 
transitions for vortex matter in type-II superconductors. 
In the stripe system the order-disorder transition occurs 
at different disorder strengths for the two different driv- 
ing directions, producing regimes of enhanced anisotropy 
in which the system depins plastically in one direction 
but elastically in the other. 

In the plastic flow regime near depinning, we observe 
a series of velocity jumps and transitions which corre- 
spond to transitions between filamentary flow states as- 
sociated with narrow band time-of-flight velocity noise 
signatures and strongly fluctuating plastic flow states ex- 
hibiting broad band noise signatures. These transitions 
are very similar to recent experimental observations in 
this class of system. 

The anisotropic transport can be enhanced by thermal 
fluctuations. Thermal disorder induces an anisotropic 
melting of the stripes, with a lower temperature stripe 
liquid in which particles can move easily along the length 
of the stripe but remain confined perpendicular to the 
stripe, and a higher temperature isotropic liquid in 
which the stripe structure is destroyed. For intermedi- 
ate quenched disorder strength, the stripe structure is 
disordered at low temperatures but can undergo a ther- 
mally induced stripe ordering into a floating stripe phase 
when the thermal fluctuations reduce the effectiveness of 
the quenched disorder. 

We expect that these results should be generic to 
any type of stripe forming system driven over quenched 
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disorder. Particular systems where the pulse measure- 
ments and transient times could be analyzed include two- 
dimensional electron gasses or the recently studied type- 
1.5 superconductors in which the vortices interact via 
competing repulsive and attractive interactions. Other 
relevant systems include stripe or labyrinth patterns in 



soft matter systems driven with electric or magnetic fields 
over a rough surface or through obstacle arrays. 

This work was carried out under the auspices of the 
NNSA of the U.S. DoE at LANL under Contract No. 
DE-AC52-06NA25396. 
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